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ABSTRACT 

Detecting heavy atoms in the inflated atmospheres of giant exoplanets that orbit close to their parent stars is a key factor for understand- 
ing their bulk composition, their evolution, and the processes that drive their expansion and interaction with the impinging stellar wind. 
Unfortunately, very few detections have been made thus far. Here, we use archive data obtained with the Cosmic Origins Spectrograph 
onboard the Hubble Space Telescope to report an absorption of ~ 6.4% ± 1.8% by neutral oxygen during the HD 189733b transit. 
Using published results from a simple hydrodynamic model of HD 189733b, and assuming a mean temperature of ~ (8 - 12) x 10 3 K 
for the upper atmosphere of the exoplanet, a mean vertical integrated OI density column of ~ 8 x 10 15 cirT 2 produces only a 3.5% 
attenuation transit. Much like the case of the hot- Jupiter HD 209458b, super-solar abundances and/or super-thermal broadening of the 
absorption lines are required to fit the deep transit drop-off observed in most far-ultraviolet lines. We also report evidence of short- 
time variability in the measured stellar flux, a variability that we analyze using time series derived from the time-tagged exposures, 
which we then compare to solar flaring activity. In that frame, we find that non-statistical uncertainties in the measured fluxes are not 
negligible, which calls for caution when reporting transit absorptions. Despite cumulative uncertainties that originate from variability 
in the stellar and sky background signals and in the instrument response, we also show a possible detection for both a transit and 
early-ingress absorption in the ion CII 133.5 nm lines. If confirmed, this would be the second exoplanet for which an early ingress 
absorption is reported. In contrast, such an early ingress signature is not detected for neutral OI. Assuming the HD 189733b magne- 
tosphere to be at the origin of the early absorption, we use the Parker model for the stellar wind and a particle-in-cell code for the 
magnetosphere to show that its orientation should be deflected ~ 10 — 30° from the planet-star line, while its nose's position should 
be at least ~ 16.7 R /; upstream of the exoplanet in order to fit the C II transit light curve. The derived stand-off distance is consistent 
with a surface magnetic field strength of ~ 5.3 Gauss for the exoplanet, and a supersonic stellar wind impinging at ~ 250km/s, with 
a temperature of 1 .2 x 10 5 K and a density ~ 6.3 x 10 6 crrr 3 at the planetary orbit, yet the fit is not unique. 

Key words, planet-star interactions- planets and satellites: magnetic fields - stars: individual: HD 189733- star: activity ultraviolet - 
sun: activity - ultraviolet: general 

1 . Introduction Two main difficulties, the relative faintness of the FUV (1 15- 

. . , „ IT , TN , . 143 nm) stellar emissions and the variability of the sources, are 

^™?cT ? ^ fa ;- ult ™ let FUV) s P ectral ™ige of at the origin of the major uncertainties on any estimation of the 
O HD 209458b, the first detected exoplanet transiting a sun-like osition and the net mass loss from a eX oplanetary atmo- 

Cn ; star, have clearly sh ow^Jhat its atmosphere uhrt and inflated here Indeed f ft 854 lanets detected until the end of 
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? P t . ga , S .° US t ° be ( |Be "' J 6 11 ? 08|) - I g be detected by the Hubble Space Telescope (HST) without re- 
the Roche lobe is a necessary but not sufficient condition for es- M expensive resources. This explains that up to now, 
b ; cape to occur since vanous effects could confine the extended om a few transiti lanets have been observed in the FIJV 
nebula. Planetary magnetic field, even of weak strength, will Furthermore> the two most aften observed starS; namel the sun . 
trap particles on closed field lines as s hown by several stud- ^ HD209458 star and the K . star HD 189m show 
les ( |Adams || 201 l H jrammell et al | 201 1|) . The impinging stellar ^ eyidence of variab0it in their flux on time scales ran 
wind can also confine the expanding g as (Bisik al o et; ing from a few h ours to the rotational period dBen-Jaffel 112007b 
Murray-Clay et al. Ett fll Stone & Praga E flflfl In addition the frj^eretriJ E^ . Flaring activity has also been reported 
planet ary magnetic field could affect ions and neutrals differ- for HD 18973 3 based on nea r i v simultaneous X-ray and FUV 
ently (j Khodachenko et al | 2012|) . The level of stellar X-ray and observations dLecavelier et al. II2012T) . A third difficulty, due to 
extreme-ultraviolet (XUV) radiation at the planet as well as time si j contamination b both the sk background signal (at some 
variations ca n dso affec t the extension of the neutral versus ion ^ ^ and ^ instrument response , adds to the complex- 
components fiio | 2011| ). All these effects cast doubts on the dis- . Qf ^ blem WMe ^ firgt difficm can bg resolyed fe 
tnbution of the upper atmospheric components assumed to pro- focusi Qn a few dose _ b and uv . bri ht starSi the si j vari . 
duce the transit absorption signatures, and thus o n the mass loss abm from both ^ source and me instrument is a rea i Wem 
rates disseminated thus lar in the literature (e.g., Lammer et al. ^ should be addressed t0 build a reliable dia „ nostic t0 extract 
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Table 1. Archival HST COS data set used in this study 



Data set name 


Exposure time (s) 


Orbital phase 


HST orbit 


IbOukq 


209 


-0.071 





IbOumq 


889 


-0.050 


1 


IbOuoq 


889 


-0.045 


1 


Ibauqq 


889 


-0.020 


2 


Ibausq 


889 


-0.015 


2 


Ibauuq 


889 


-0.009 


2 


IbavOq 


889 


+0.010 


3 


lbav2q 


889 


+0.015 


3 


lbav4q 


889 


+0.021 


3 


Ibavdq 


889 


+0.040 


4 


Ibavmq 


889 


+0.045 


4 


Ibavpq 


889 


+0.051 


4 



Notes. The twelve individual exposures were obtained with the G130M grating in the timetag mode (pr ogram HST-GO-1 16 73). Central transit time 
was defined by propagating from zero phase on HJD 2453988.80339 and a period of 2.21857312 days ( Triaud et al. 2009). The short exposure #0 
has not been used in this work. 



an accurate description of the upper atmosphere and of the in- 
teraction region between the exoplanet and the impinging wind 
from its parent star. 

Before one can conclude on the relevance of the key pro- 
cesses that shape the upper atmosphere of hot-Jupiter exoplan- 
ets, it is important to know their basic composition. The FUV : 
regime is suitable for detecting atoms that show resonance lines ° 
in that spectral window. For the exoplanet HD 209458b, HI, 1 
OI, and CII were detected in the FU V, while doubts remain 
on Si III dBallester & Ben-Jaffel 1 12011 . For b oth HP 189733b 
and 55 Cnc b , only H I has been detected so far dLecavelier et al. I 
I2TMI20T1 Ehrenreich et al. 1 I20T1 . Other ato mic species ob- 
serve d to date on HD 189733b are neutral Na I dHuitson et al. I 
120 121) and K I llPont et al. \\20H) . but these species are detected 
at the base of the thermosphere and below (below ~5000km al- 
titudes) since they are readily ionized by stellar UV light. 

Here, we use archival data obtained in 2009 with the Cosmic 
Origins Spectrograph (COS) onboard HST to study the charac- 
teristics of the HD 189733 system in more details, and address 
the problem of signal variability and its effect on determinating 
atmospheric abundances and properties of the star-planet inter- 
action regions of these distant objects. To analyze the observa- 
tional results, we estimate the properties of the stellar wind at 
the planet distance and adapt an existing model to simulate the 
wind's interaction with a magnetized hot- Jupiter, tailored to the 
specific case of the HD 189733 system. 



2. Observations and data analysis 

The medium-resolution G130M grating was used and kept at a 
single wavelength setting (1291 A) throughout the observations, 
made on Sep 1 6 ( 1 8 : 3 1 UT) - 1 7 (0 1 : 1 1 UT) 2009, to best deter- 
mine relative changes in the stellar signal caused by the plane- 
tary transit (see Table 1). The co-added spectrum of the star is 
shown in Figure 1 . The table lists the exposures, orbital phases, 
and the HST orbit during which each exposure was obtained. 
Here, we use the time-tagging of the COS data to derive a time 
series of the system signal ever y 100 s, following the procedure 
described in iBen-Jaffell (12007). This allows us to better deter- 
mine the true transit absorption signature and obtain a good track 
of the signal variability over selected periods. Furthermore, the 
time series should reveal any sudden stellar variation that may 
occur during short periods. We find that although the observa- 
tions did not sample the middle of the transit, they did sam- 




Fig. 1. HD 189733full spectrum obtained by HST/COS G130M. 
The geocoronal O I and H I emissions have not been subtracted. 



pie a good portion of the transit. Indeed, the exoplanet's orbital 
phases corresponding to HST orbits 2 and 3 are ins ide the transit 
windo w observed for HD 189733b in the visible dTriaud et al. I 
120091) . 

The derivation of the in-transit/out-of-transit flux ratios, and 
thus absorption by the planet disk and atmosphere, is sensitive to 
the choice of the out-of-transit reference exposures. Our choice 
here was based on the behavior of the stellar O I signal. After 
several trials, we opted to merge the last three exposures for 
the out-of-transit reference (see Table 1). This selection was also 
tested by using it to extract the transit signature from the mod- 
erately bright stellar lines Si IV 139.3 and 140.3 nm. Their in- 
tegrated flux over the full line profiles gives an obscuration of 
2.6±2.5% (statistic al), which agrees w ith the nominal disk ob- 
scuration of 2.5% dBouchv et al.ll2005l) . 

For all lines, a weak stellar continuum has been subtracted, 
and photon noise errors were propagated. Transit absorptions for 
the most relevant species are listed in Table 2. In the follow- 
ing, we focus on the four species, O I, C II, Si III, and Si IV for 
which the stellar signal is strong enough to allow a reliable study 
of the absorption during the planetary transit and/or to address 
the unavoidable variability due to all sources within the field 
of view and from the instrument response. Among the selected 
lines, only O I has a background emission from the geocorona. In 
addition, the Si III 120.6 nm emission is contaminated from the 
extended blue wing of the stellar Lyman-o' line. This is one of 
the reasons why we kept the Si III line to evaluate the time vari- 
ability, but did not use it for any diagnostic of the planetary ab- 
sorption. Light curves were first derived and studied statistically 
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Wavelength (nm) 

Fig. 2. HST/COS G130M OI 130.2 nm . Solid line: combined 
HD 189733 and background geocoronal lines obtained on egress 
(dataset lbavdq). The thin feature is the star emission and the 
broad one is the Earth line emission filling the 2.5 arsec aperture. 
Dotted line: independent spectrum of the geocoronal OI emis- 
sion line that was scaled to subtract the sky background contam- 
ination from the stellar observation. 



to evaluate the confidence we may attach to simple representa- 
tions of the time series before we compared in- to out-of-transit 
line profiles of the key lines. The crucial problem of the time 
variability of the observed flux, particularly that related to the 
patchiness on the stellar disk is discussed in section 3. 

2.1. Oxygen triplet 130.4nm 

The 100-s time-tag sampling of the exposures does not show 
any steep change in the total flux for the oxygen emission in any 
of the triplet lines beyond statistical fluctuations as a function of 
time. Significant geocoronal background emission is seen partic- 
ularly in the first exposure of the HST orbit, and we successfully 
removed this background with a method that uses a deep-sky ex- 
posure taken by the COS team in the same grating wavelength 
setting (http://www.stsci.edu/hst/cos/calibration/airglow.html). 
First, we notice that the HD 189733 O I lines appear as thin and 
bright features over the broader gecoronal lines (Fig. 2). The ad- 
ditional error from the imprecision on the geocoronal subtraction 
represents few to 30% of the total error depending on the emis- 
sion level of the contamination. Errors shown here do include 
that additional uncertainty. It is important here to stress that for 
the broader stellar Lyman-ff emission, the line width is compa- 
rable to the width of the geocoronal line that fills the science 
aperture, as seen by COS, which prevents sepa rating the stellar 
line from the background (Linsk vet al. 1 120 10h . The technique 
used here offers a unique opportunity to extend the use of COS 
to derive O I abundances from transiting exoplanets for which 
the O I lines are thin and bright enough to allow a spectral sepa- 
ration from the extended geocoronal line. 

The OI 130.4nm feature is a triplet with lines at 130.2, 
130.49, and 130.6 nm. However, the 130.49 nm line is contami- 
nated on the blue wing by the stellar Sill 1 30.44 nm line, which 
hinders subtracting the geocoronal O I line. Figure 3 shows the 
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Fig. 3. HD 189733b transit light curves obtained for selected 
spectral lines versus time (or orbital phase): OI (130.4nm) 
triplet, CII (133.5 nm) doublet, Si III (120.6 nm), and Si IV 
(139.8 nm) lines. Usually, those lines ar e used as a diagnos - 
tic of solar flux variability in the FUV dBrekke et al. I fl996). 
A synthetic light curve corresponding to an occulting disk of 
~ 1.7 R p is shown with C II and O I data points. Dark error bars 
are statistical while the gray bars include the non-statistical scat- 
ter of the stellar signal within the given HST orbit. 



transit flux ratios derived by merging the OI triplet lines, inte- 
grating each line out to ~ ±40 km/s, since farther out there is sig- 
nificant error from the subtracted geocoronal background. A chi- 
square analysis of the OI transit light curve shows that its repre- 
sentation by a linear trend, which could be related to a patchiness 
effect or a temporal variability in the stellar emission or a combi- 
nation of these, is rejected with a confidence level slightly higher 
than 99% (10 degrees of freedom). This result was obtained us- 
ing the PAN package that implements the Levenberg-Marquardt 
technique and a bootstra p Monte Carlo simulation of the statis- 
tical noise propagation (Dimeo l l2005h . In contrast, assuming a 
synthetic light curve model (occulting disk with the radius and 
the time center as free parameters) gives ax 2 ~ 12(10 degrees 
of freedom), a good indication that the O I dataset is consistent 
with a transit by an occulting disk of ~ 1 .7 R p radius. This result 
is also confirmed by the comparison of the in-transit (orbit 2) 
to out-of-transit (orbit 4) line profiles for the OI 130.6nm line 
(Figure 4a). The absorption derived for the 130.2 and 130.6nm 
lines combined is ~ 8.1 ± 2.2%, while that for the three lines 
is ~ 6.4 + 1.8%. We adopted the ~ 6.4 + 1.8% absorption as 
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the proper measurement. Despite the noise level and the limited 
accuracy on COS/G130M wavelengths, a transit absorption line 
width of ~ 50± lOkm/s is estimated for both the 130.2 and 130.6 
nm lines after correcting for the line spread function (LSF). 

It is important to stress that the reported uncertainties in- 
clude photon noise but no stellar/detector variability (see Table 
2 for more details). The size ~ 1.7 R p of the OI opaque region 
falls inside the 2.84 R p Roche lobe radius (which would pro- 
duce 19% obscuration), yet this does not mean that the Roche 
lobe is not filled. Indeed, the intense X-ray and EUV flux from 
the young and active K1V star may easily ionize the gas in- 
side t he external layer between 1.7R ; , and 2.84R p that faces the 
star (IGuo 11201 II) . More details on the ionized components are 
provided in the following section when we analyze the singly- 
ionized carbon lines. Interestingly, the low-resolution detection 
of extended H I on this planet has shown a confined ~ 5.0 + 0.8% 
absorption averaged oy er three transits observed in 2007-2008 
(Lecavelie r et al. 1 12010). Recent medium-resolution HST/STIS 
H I Lyman-a transit data seems to indicate a variation with val- 
ues of 2.9 ± 1 .4% in 2 010 and 5.0 + 1.3% in 201 1 for the line- 
integrated absorption dLecavelier et al. 1 120121) . Our transit ab- 
sorption obtained for neutral oxygen is thus consistent with most 
neutral hydrogen detections, despite the reported variations. The 
time variation problem for the COS data set is discussed in more 
detail in Sect. 3. 

2.2. Carbon doublet 133.5 nm 

This stellar emission is the brightest FUV emission after HI 
Lyman-a (and about as intense as the HI Lyman-a observed 
from HD 209458), and thus provides relatively high signal-to- 
noise (S/N) data. After subtracting a weak stellar continuum 
level, the ratios of the stellar flux per exposure versus the nom- 
inal out-of-transit data were derived by merging the CII 133.4 
and 133.6 nm lines, integrating each line out to ~ ±50 km/s from 
line center. These ratios are shown in Figure 3. A chi-square 
analysis of the C II time series leads to a x 2 ~ 40.0 when us- 
ing a linear trend (such as related to patchiness and/or temporal 
variation; 10 degrees of freedom), while for a representation by 
two synthetic light curves, we obtain ax 2 ~ 28.8 (3 free param- 
eters: depth of planet's transit, depth and relative phase position 
of an extra absorption such as a magnetosphere; 9 degrees of 
freedom). The better confidence obtained for the two-light-curve 
models may not be real because the extra freedom allowed by the 
model may enhance the probability that a false-positive result 
would be obtained. Nevertheless, the statistical analysis supports 
the conclusion that both the linear trend and the two-light-curve 
models poorly represent the data set. Obviously, a new approach 
that goes beyond the pure statistical analysis is required to im- 
prove our understanding of the C II timeseries. 

In a first step, we compared the CII 133.5 nm line profiles 
respectively from the in-transit and ingress orbital phases to the 
out-of-transit line profile. As shown in Figure 4b, when com- 
pared to orbit 4, the orbit 2 (in-transit) line profile has a clear 
absorption that is higher than the statistical noise. Based on this 
CII line profile diagnostic, a transit absorption of ~ 7.0 ±2.1% 
is derived. Despite the noise level and the limited accuracy on 
COS wavelengths, we can also estimate that the spectral width 
of the transit absorption line is ~ 60 km/s after correcting for the 
instrument's LSF. Similar to the transit absorption obtained for 
O I, the C II absorption would correspond to an opaque occult- 
ing disk of ~ 1 .7 R p that is smaller than the Roche lobe. The C II 
fluxes also show a significant absorption, at the ~ 10% level, ear- 
lier in phase (during orbit 1) than expected from ions confined to 
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Fig. 4. Comparison of in-transit (orbit 2 : dotted) to out-of-transit 
(orbit 4: solid) line profiles. We binned spectra by 4 (~ 7.5 km/s), 
approximately corresponding to the nominal spectral accuracy 
of the dispersion of the COS/G130M grating. A redshift of 
~ 2.6 km/s of the star was corrected for (Bouchv et al. 1120051) . It 
is important to stress that the COS wavelength calibration does 
not allow one to know the relativ e position of the line s with an 
accuracy better than ~ 7.5 km/s (Oliv eira et al. Il201~0h . For the 
available low S/N COS dataset, this inaccuracy makes the in- 
terpretation of the transit absorption line profile very uncertain, 
particularly the veracity of any blueshift or redshift that may 
appear between in-transit and out-of-transit line profiles, (a) OI 
130.6 nm line. (b)CII 133.5 nm line, (c) Si IV 139.3 nm line. 
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the Roche lobe (e.g., Fig 3). The early-ingress absorption is fur- 
ther confirmed by the line profile diagnostic (e.g., Figure 5). This 
result suggests properties in the C II ions that are strikingly dif- 
ferent in the O I neutrals in the extended upper atmosphere and 
beyond, where plasma interactions are at play. It is unlikely that 
the apparent extra absorption is caused by intrinsic stellar varia- 
tion since no other stellar emission shows this temporal behav- 
ior at high S/N. The temporal variation problem is discussed in 
more detail in the following section to check the veracity of the 
observed event. If confirmed, this would be the second target af- 
ter WASP- 12b for w hich an early-ingress absorption is detected 
dFossatiet al. 112010b . 

In the following, we therefore address the time variabil- 
ity problem of the measured fluxes, focusing on the different 
sources that may produce it. To distinguish the different effects, 
we propose a multi-species (spectral) analysis that should clar- 
ify most uncertainties and confirm our diagnostic about the O I 
detection and the potential C II detection. 



3. Temporal variation: a combined statistical and 
multi-spectral approach 

Addressing the temporal variation problem during a planetary 
transit is very complex because the observed variability could 
originate from the full stellar disk, the stellar latitudes sampled 
during transit, the details of the planet-obscuring area, and the 
instrument. Each of these sources of variability represents a chal- 
lenge in itself, particularly if the data set is limited in time, which 
offers no opportunity to check the repeatability of any event. In 
the following, we discuss the three sources of variability as they 
may appear in the COS dataset, focusing on emissions of the four 
species O I, C II, Si III and Si IV. Our goal is to obtain a higher 
confidence in the OI and CII detections, particularly by checking 
if a stellar variation might mimic a transit-like light curve in the 
current COS dataset. 

3.1. Stellar intrinsic activity 

Stellar activity is a crucial problem in characterizating extra- 
solar planets and their evolution. The activity of HD 189733 
in the X-ray and EUV has been reported from XMM-Newton 
X-ray spectra and from modeling of the EUV emission that 
was based on prev i ous s tudies of stellar coronal models by 
ISanz-Forcada et al71 (1201 lb . The luminosities are log Lx = 28.18 
and logL EUV = 28.48 in the X-ray (0.5-10nm) and EUV (10- 
92 nm). These values are almost identical to those of the ex- 
tensively studied star Epsilon Eridanis (logLx = 28.20 and 
logL EUV = 28.44). Epsilon Eridanis and HD 189733 are both 
early-K main-sequence stars (K1-K2 and K0-K1) of relatively 
high acti vity and similar est i mated ages (respectively 1 . 1 and 
1.2 Gyr; ISanz-Forcada et afTI (1201 lh ). X-ray measurements of 
HD 189733 have also been made with the SWIFT satellite span- 
ning about 30 hours in 201 1, during wh ich time enhanced activ - 
ity was seen, including a bright flare dLecavelier et al. Il2012h . 
The average of the combined X-ray and EUV fluxes derived 
from the SWIFT data was 7.1xl0 28 ergs/s. This is somewhat 
higher than the total X-ray and EUV 4.5 xlO 28 ergs/s fluxes 
previously reported (ISanz-Forcada et a"Tl 1201 lh . The chromo- 
spheric Call H & K line emission of HD 1 89733 has also 
been measured at log(/?Hic)=-4.501 with Keck dKnutson et al. I 
l2010h . The activity of HD 189733 is higher than that of 
the Sun. The me an Sun has logLx - 27.35 erg/s and 
log(/? HK )=-4.905 (Mamaiek &Hillenbrandl 12008: Mamajek 



120121: lJudge et al. I l2003h . For comparison, HD 209458 may 
be of medium, sun-like activity, although there are significant 
uncertainti es (Balle ster & Ben-Jaffel 1120131 : ISanz-Forcada et al. I 
1201 lUKnutson et al. 11201 0h ~ 

Available FUV observations of HD 189733 are very lim- 
ited in time and offer no full reference on the activity level 
versus short and long time scales from existing data. What re- 
mains are the solar activity databases as a reference for the rel- 
ative variation of the key FUV resonance emission lines. Our 
goal is to derive a reliable diagnostic on the activity level that 
occurs during the observing period reported in this study, us- 
ing the solar flaring activity as a reference. For this purpose, 
we revisited Solar Radiation and Climate Experiment (SORCE) 
FUV archive observations of the Sun obtained since 2003 along 
with a catalog of flares observed by the Solar EUV Experiment 
on the Thermosph ere Ionosphere Me sosphere Energetics and 
Dynamics mission ( Wood s"et al. 1 12006). 

Without any loss of generality, we may classify the stel- 
lar activity into long-term (longer than spin period), the mid- 
term (close to a spin period), and short-term (a day or less). 
Furthermore, in addition to the solar periodicity related to both 
the 27-day rotation and the solar cycle, flaring activity may 
occur at a level that c an be, for strong events, comparable to 
the long-term activity dBrekke et al. II 19961 IWoods et al~ll2003l) . 
The variation in Si III, Si IV, HI, O I and C II FUV lines have 
been reported for s unspot activity ([Woods & Rottmanl [2002; 
IWoods et al.ll2006t ISnow et al. Il2010l) . However, for flare ac- 
tivity, only very strong or moderate events have been pub- 
lished that we are aware of (Bre kke et al. 1119961 IWoods et al. I 
2003). To obtain several levels of activity, namely low, mod- 
erate, strong, and extreme levels of flare activity, we analyzed 
SORCE data obtained between 2003 and 2007, a period for 
which a flare activity catalog exists. Solar spectra obtained with 
the Solar Stellar Irradiance Experiment (SOLSTICE) onboard 
SORCE on a daily basis have a low 1 nm resolution, yet they 
offer the opportunity to study most of our resonance lines si- 
multaneously. After subtracting a weak solar continuum, we de- 
rived a timeseries of the irradiance of each line during the se- 
lected period. Using the list of flare events provided in the LASP 
database (lasp.colorado.edu/see/see-flare-catalog.html), we de- 
rived for each line an approximate estimate of the activity level 
that is calculated as a standard deviation around each event de- 
scribed in the LASP catalog (see Table 3). It is interesting to find 
that the variation levels are comparable for all lines for the low- 
activity sector, which may include variation of the plages and 
enhanced networks. Results differ significantly for stronger flare 
events. For extreme level flares (X17), the stellar variation is 
strong but remains comparable to the sunspot activity (see Table 
3). 

Here, we emphasize that the solar activity derived was used 
as a reference level for comparison with the HD 189733 activity 
for the same resonance lines. We did not use these numbers to 
derive any intrinsic activity level of the HD 189733 FUV line 
emissions. We instead preferred to rely on the HST/COS data 
thus far obtained in the time-tag mode to derive a time series 
that was used to estimate stellar/instrument temporal variations 
that were then added quadratically to the statistical noise. In this 
way, the shown light curves have more realistic error bars. 

3.2. FUV HD 189733 flux variability: assessment 

To track stellar activity during the transit observations of 
HD 189733b, we generated a time series of observed spectra ev- 
ery 100s. The time series has several gaps, yet it is useful for 
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Table 2. Absorption depth of stellar emission during the HD 189733b transit 

Atoms Wavelength (nm) Flux drop-off " (%) Photon-noise error (%) Total error bar (%) 



OI 130.21 + 130.60 8.1 2.2 3.4 

OI 130.4 triplet 6.4 1.8 3.4 

CII 133.45 + 133.57 7.0 1.2 2.1 

Si III 120.6 6.8 2.0 4.7 

SilV 139.3 +140.3 2.6 2.5 5.6 



Notes. Flux drop-off during transit of HD 189733 derived from HST/COS G130M medium-resolution observations. Total error bars include both 
photon noise and stellar/instrument variability, (a) The signal drop-off is derived from the data at planet orbital phase of ~ 0.985 sampled during, 
transit in HST orbit 2. 
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same as in Fig. 4. (a) CII 133.5 nm line, (b) Si IV 139.3 nm line. 



deriving key information versus the HST orbit for the spectral 
line selected. To this end, we first extracted the time evolution 
of the merging of the OI 130.4nm triplet lines, the merging of 
the CII 133.5 nm doublet, the merging of the Si IV 139.8 nm 
doublet, and the Si III 120.6 nm line (Fig. 3). In a second step, 
for each HST orbit and for each species, we gathered data every 
~ 900 s and derived the corresponding scatter observed between 
the resulting spectra during that orbit. The signal scatter was es- 



timated as a standard deviation from the mean flux of the desired 
resonance line at the selected HST orbit (or transit phase angle). 
Finally, the photon noise level attached to the mean flux obtained 
during the HST orbit was subtracted quadratically from that spe- 
cific standard deviation between the ~ 900s sub-spectra. The 
resulting variability in the signal level could be of stellar, plane- 
tary, or instrumental origin. To facilitate the interpretation of this 
variability, we plotted that extra variability versus the HST orbit 
for the four resonance lines selected in Figure 6. In the same 
figure, we also show the Sun's medium flaring activity previ- 
ously derived in Table 3 as a reference. We note that a variability 
zero means that the signal is dominated by the photon noise. 
Representative photon-noise error bars are shown in Figure 6 at- 
tached to the solar value corresponding to each line. 

First, we notice that most temporal variations (beyond pho- 
ton noise) recorded during the total duration of the observations 
(~ 5 hours) are weaker than the Sun's medium flaring activ- 
ity. Only for orbit 3, the extra variability recorded for the Si IV 
139.8 nm lines is comparable to the Sun's medium flaring ac- 
tivity. For that orbit, the 900 s data scatter seems to indicate a 
relatively strong temporal variability that is clearly visible for 
the Si III 120.6 nm and Si IV 139.8 nm lines. However, for the 
O I & C II lines, the derived scatter of the ~ 900 s signal seems 
smaller for orbit 3 than for the Sun medium flaring activity, yet 
it is higher than the photon noise level. This specific behavior 
of the orbit 3 signal will be very useful in the following for dis- 
cussing the different sources of the observed fluctuations. The 
similitude between signal scatter levels observed for the O I and 
C II lines shows that the origin is probably not the residuals from 
the correction of the sky background temporal variation because 
the C II lines are not affected by that contamination. For a total 
exposure time of ~ 45 minutes, we are also at a loss to explain 
the high variance by the response of the COS detector at the 
Si IV spectral position. In addition, as we discuss below, for the 
O I and C II lines, the crossing of a bright local region by the 
planetary obscuring area does not seem to show strong fluctua- 
tions during the 45 minutes of observation during orbit 2. This is 
demonstrated in figure 4, where we see for orbit 2 (during tran- 
sit) comparable or lower fluctuation than in orbit 4 (out-of-transit 
reference orbit) for the OI and CII lines. For all these reasons, 
we can conclude that the signal variations observed for OI, CII, 
and SI IV on the orbit-to-orbit time scale are most likely of stel- 
lar and/or planetary origin (e.g., Fig. 6). The key question is how 
to distinguish between the two fluctuation sources on the basis 
of the limited COS dataset? 

To begin answering this question, we stress that chromo- 
spheric OI and CII emissions are not restricted to the active 
regions only (the so-called patches). There are background- 
quiet regions that do emit in the FUV lines that originate from 
the chromosphere but are not that dependent on solar cycle 



6 



Ben-Jaffel & Ballester: Detection of oxygen in exoplanet HD 189733b 



< 



15 



S 10 
o 

01 



> 



Si III 



* 



▲ 

o 
¥ 



HD 189733: FUV Activity 

-1 | i i i i | i i i i | 

01 Cll 

i I 



Si IV _ 



* 



V o 



1200 1250 1300 1350 
Wavelength (A) 



I 400 



Fig. 6. HD 189733 activity at key species emission lines versus 
HST orbits during the 5-hour duration of the HST/COS G130M 
program used in this study. The solar flare activity is taken from 
Table 3. Open circles denote HST orbit 1, inverted open triangles 
orbit 2, filled triangles orbit 3, diamonds orbit 4, and stars the 
solar medium flare activity. The corresponding statistical noise 
levels are shown as gray error bars. 



dWoods et al. ll200(AlWo"rden et al. I200lh . Since we do not know 
the quiet versus patchy distribution on our particular star, we 
cannot estimate the exact effect a particular bright patch cross- 
ing would have. However, we do know that the patchiness on 
the star is more pronounced for these emission lines originat- 
ing from the hotter regions. Taking the Sun as a reference, we 
remark that the Si IV 139.7 nm doublet originates from a hotter 
layer at log(T)~ 4.75 tha n O I 130.4nm (log (T)~ 3.85) and C II 
133.5 nm (log(T)~ 4.1) (I Woods et al. 1 2000). First, we need to 
verify whether what is known for the Sun is consistent with the 
HD 189733 COS observations. Recalling Fig. 3, we see that the 
O I and C II lines have much less scatter within orbits than the 
Si IV lines if the entire intra-orbit variation in orbit 3 and 4 were 
generated from patchiness on the star. This result is quite consis- 
tent with FUV observations of the Sun where the contrast in the 
plages and enhanced network regions is stronger for the Si IV 
139.3 nm and Si III 120.6 n m lines than for the OI 130.4nm 
and C II 133.5 nm emissions (IWorden et al. l200lh . This result is 
very important because it allows us to use the comparison be- 
tween O I & C II from one side and Si IV lines from the other 
side to now investigate patchiness/transit effects from one orbit 
to another. 



3.3. Temporal variation: final diagnostic 

For the OI data, our chi-square analysis rejected a linear trend 
with a very high confidence. In addition, the same statistical 
analysis showed that a transit light curve (with the same degrees 
of freedom as the linear trend) was a good representation of the 
data. Furthermore, the comparison between in- to out-of-transit 
line profiles for the O I lines clearly shows the expected absorp- 
tion. Finally, the observed flux drop-off in the O I signal that ap- 
pears unchanged between orbit 2 and 3 is unlikely to be produced 
by the planet crossing a concentration of active regions because 
we see no similar effect in the Si IV lines which are more sen- 
sitive to these regions. Indeed, we emphasize that for the Si IV 



data set, the orbit 2 average signal does show the expected tran- 
sit absorption, while orbit 3 shows the opposite effect from that 
observed in the O I signal. We can therefore conclude that the 
O I transit is detected with a high confidence. 

For the C II dataset, assuming a linear trend is a good repre- 
sentation of the time variation of the signal (despite the poor con- 
fidence derived from the statistical analysis), one would expect 
a more pronounced effect for Si IV emission because the latter 
originates from hotter regions. As shown in Fig. 3, the orbit-to- 
orbit variation in the Si IV averaged-signal is much weaker than 
the one observed for C II. In addition, the flux averaged in or- 
bit 2 shows a variation that is exactly what is expected from a 
transit for the Si IV line. For these reasons, we dismiss the idea 
that any spatial concentration of active regions on the stellar disk 
may have produced the C II light curve. The fact that the scatter 
is larger for Si IV lines during orbits 2 and 3 is probably related 
to bright-limb-crossing effects. Indeed, the expected limb bright- 
ening for Si IV (a center-to-limb factor (CTL) ~ 3.82 increase in 
the intensity for the Sun) is much higher than for t he O I and C II 
lines (CTL~ 1.21 for CII, CTL~ 1.19 for OI) (IWorden et al. I 
l200lb . For Si III the bright-limb -crossing effect is also signifi- 
cant with CTL~ 1 .74. In addition, the timing of the two orbits (2 
and 3) is close to the planetary disk crossing of the stellar bright 
limb. Therefore, the scatter thus far observed (high, low, high 
flux) for Si IV during orbit 3 is probably related to the bright- 
limb-crossing with the middle point being exactly the crossing 
time on egress. The bright-limb-crossing effect also seems to be 
happening in orbit 2. 

The averaged signal from orbit 3 shows for C II a specific 
behavior that is observed in all lines except O I, and that is not 
easy to explain with simple changes in the planet absorption or 
by the stellar patchiness. The origin of this enhancement could 
be a flare event, a scenario that cannot be assessed from the lim- 
ited dataset however. Despite the uncertainty on the origin of the 
orbit 3 signal, the C II line profile information from orbits 2 and 
4 and the inconsistency of the patchiness assumption discussed 
above (see Figure 5a) both support a tentative detection of the 
C II transit. In addition, because we do not see any similar behav- 
ior in the Si IV lines that would result from a patchiness effect, 
we now have more confidence in the early-ingress absorption 
possibility for C II, particularly in view of the line profile com- 
parison in Figure 5. Nevertheless, in all cases, we must stress 
that results for the C II data set require confirmation with future 
observations. The tentative detection of an early-ingress ion ab- 
sorption on this planet has motivated us to model the interaction 
of the planet with the magnetized stellar wind. This modeling 
effort is a valuable tool for planning future observations of the 
system. 

4. Comparison with model simulations 

The O I lines show an ~ 7 % level of absorption during transit. 
This absorption can be caused by an opaque disk with a size of 
~ 1.71^, that is smaller than the Roche lobe size (~ 2.84 R^ 
along the planet-star line and ~ 2.06 R p across) of the exo- 
planet. To properly constrain the atmospheric composition, so- 
phisticated models that include photochemistry and radial trans- 
port in the atmosphere of HD 189733b are needed. In addition, 
a careful analysis of the transit absorption requires a parametric 
study that is beyond the scope of this paper. However, using sim- 
plified modeling that accounts for the hydrodynamical transport 
and ionization of the main species in the extended atmosphere 
is enough for our purpose to achieve a better insight into the O I 
density column required to fit the observed transit drop-off of 



7 



Ben-Jaffel & Ballester: Detection of oxygen in exoplanet HD 189733b 
Table 3. Solar flare activity in the FUV 



Species (feature nm) 


Activity level 


Nature 


Variability level" (%) 


CII (133.5) 


(L, M, H, E) 


flare 


(3.7, 7.6, 12.5, 30.0) 


01(130.5) 


(L, M, H, E) 


flare 


(2.7, 5.0, 7.3, 16.0) 


Si 111(120.6) 


(L, M, H, E) 


flare 


(4.0, 9.0, 12.8, 27.6) 


Si IV (139.8) 


(L, M, H, E) 


flare 


(4.3, 8.6, 14.3, 33.4) 


CIV (154.0) 


(L, M, H, E) 


flare 


(3.9,5.0, 11.7,31.7) 


CII (133.5) 


(27D, 11Y) 


sunspot 


(22, 57) 


01(130.5) 


(27D, 11Y) 


sunspot 


(15, 29) 


Si III (120.6) 


(27D, 11Y) 


sunspot 


(27, 73) 


Si IV (139.3) 


(27D, 11Y) 


sunspot 


(22, 60) 



Notes. Solar observations were obtained with the Solstice instrument onboard SORCE using the 1 nm low-resolution spectrometer. The flare 
activity level is based on a list of flaring events derived by the Colorado/LASP team (lasp.colorado.e du/see/see-flar e-catalog.html). (L, M, H, E) 
are for low- (class Ml), me dium- (class M5), high- (X I), and extreme- (X17) activity flare activity (Hudson 2011). (27D, 11Y) are the 27-days 
spin and solar cycle periods (Woods & Rottman 120021) . (a) The variability level is estimated as a standard deviation around the selected flare event. 
For the 11Y period, the variability is the ratio of the mean solar irradiance value when the solar cycle is at its maximum (1992) to its minimum 
(1996). 



~ 7%. Recent photochemistry and hydrodynamic models show 
that the base of the planetary wind is close to t he pressure level 
of ~ 0.1 - l^ubar at whi ch H2 is dissociated dKoskinen et al. I 
I20K1 iMoses eTalllloTl . In the following, we assume that the 
atmosphere is spherically symmetric and that OI is dragged out 
to high altitudes by H I dGarcfa Munoz I l2007t iKoskinen et al. I 
2012a). This allows us to build the atomic oxygen distribution 
versus altitude from the H I distribution embedded between the 
bottom of the planetary wind region and the edge of the Roche 
lobe, assuming a hydrostatic regime for the atmosphere dGuo I 
1201 ll) . For reference, the HI and temperature distributions used 
here are based on a multi-fluid model of atomic-proton mi xture, 
yet it does not include photochemistry nor heavy atoms dGuo I 
2011). In a first step, we assume that only thermal broaden- 
ing affects the atmospheric absorpti on. Starting from the solar 
abundance distribution provided by iGuo I (1201 ll) . we can esti- 
mate an ~ 3.5% attenuation during transit that corresponds to 
OI density column of ~ 8.3 x 10 15 cm 2 , and a volume den- 
sity noi ~ 6.8 x 10 7 cm~ 3 at the base of the model near r=R p . 
Scaling the initial atmospheric model to a 40x solar-abundance 
level produces an attenuation of ~ 5.3% during transit. The cor- 
responding OI density column is ~ 3.3 x 10 17 cm 2 produced 
with a volume density hqi ~ 2.7 x 10 9 cm~ 3 at r = R p . These 
high densities are probably difficult to justify because that would 
increase the atmospheric mean mass, descrease the correspond- 
ing scale height, and consequently would require unrealistically 
high er stellar EUV input to m aintain the extended atmosphere 
(IKoskinen et al. l[2010L[2012bl) . 

In a second step, we include super-thermal line broadening 
that is produced by a layer of hot O I atoms c onfined on top of 
the atmos phere (iBen-Jaffel & Hosseinill2010l) . As described in 
detail in IBen-Jaffel & Ho sseini (120101) . both the bottom posi- 
tion and the effective temperature of the layer are free parame- 
ters. Our intent is not to conduct a parametric study, but rather 
to investigate how super-thermal populations may modify the re- 
sult obtained above for only thermal brodening. Assuming solar 
abundances, a hot OI layer above ~ 1.1 R p with an effective 
temperature of Toi ~ 8.4 x 10 4 K produces a transit absorption 
of ~ 4.6% that is ~ lcr below the nominal value. Increasing the 
solar abundance by x5 solar for the same hot O I layer produces 
the observed ~ 6.4% transit drop-off. If we place the bottom of 
the hot layer near ~ 1.5 R ;; , the same 5x super-solar model with 
the same effective temperature of ~ 8.4 x 10 4 K leads to a transit 



attenuation of ~ 5.2%. Other solutions exist and may require a 
full sensitivity study. 

It is important to stress that our results are not senstitive to 
the OI line selected for the atmospheric analysis. In agreement 
with recent results on HD 209458b, super-solar abundances with 
thermal broadening and super-thermal b roadening probably af- 
fect the FUV transit of hot e xoplanets (IBen-Jaffel & Hosseini I 
l2010t IKoskinen et al. 1 12012b). Sophisticated models that com- 
bine photochemistry and hydrodynamic transport in the atmo- 
sphere of HD 189733b are needed in the future to test whether 
solar or super-solar abundances in the exoplanet's bulk atmo- 
sphere may produce enough thermal a nd super-thermal OI 
atoms to fit the den sities derived here dGarcfa Munoz 1 120071 
IKoskinen et al. I l2012al) . In addition to a sensitivity study that 
includes most important line brodenings, future observation of 
high-resolution line profiles of key FUV lines must be obtained 
with high S/N to separate the different effects. 

In contrast to the neutral O I component for which a classical 
transit light curve is obtained, the peculiar signature detected in 
the early ingress of the C II lines seems to indicate that the gas 
topology around the very XUV-hot Jupiter HD 189733b cannot 
be sketched by the Roche lobe volume simply filled by the mag- 
netized plasma. Future observations are much needed to confirm 
the veracity of the C II light curve, yet we can already explore 
whether such a feature could be predicted by model simulations 
using realistic parameters for the HD 189733 system. 

4.1. Model simulation of HD 189733b magnetosphere 

An early-ingress absorption feature was alread y reported in 
the n e ar-UV for the very-ho t-Jupiter WASP- 12b dFossati et al. I 
120101: lHaswell et al. I 120121) . The first explanation proposed 
thus far is that the extra absorption is produced by the 
denser and hotter gas of the bow- shock upstream of the ex- 
planet magnetosphe re (Fossa tTet al. 11201 Ol IVidotto et al. fe OlO: 



Khodachenko etaT] 120121) . To study the stellar wind interac- 
tion with the exoplanet, a few numerical si mulations have been 
made using MHD and hybrid 3D codes (jP reusse et al. 2007; 
iJohansson et al. I l201lt ICohen et al.l l201ll) . While large-scale 
simulations that include the star-planet system are very useful 
for obtaining a global view of the system, they do not give that 
much detail on the planetary magnetosphere itself or on its opac- 
ity imprint on the transit light curves. In addition, the intrinsic ro- 
tation of the magnetized planet and the formation of the magne- 
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tospheric current sheet are not properly inco rporated in those 3D 
nume rical simulations (for more details, see Khodach enko et al~l 
l20T2h , Therefore, up to now, only empirical models have been 
used to derive the bow-shock positi on and predict its opac- 
ity effect on the transit light curves dKhodachenko et al. 11201 2t 
iLlama et al. ll201lllTrammell et al. ll201ll )~ 

Another explanation of the early-ingress absorption relies on 
gas-dynamics simulations of the stellar wind interaction with the 
exoplanet's atmosphere where no magnetic fields are included. 
These simulations predict a large-scale asymmetric gas distribu- 
tion around an unmagnetized exoplanet tha t is supposed to pro - 
duce a distorted light curve during transit dBisikalo et al. 112013b . 
While the stellar wind and planetary magnetic fields are difficult 
to estimate in general, the gas-dynamics models may produce 
very useful scenarios of stellar interaction for very weakly mag- 
netized exoplanets, particularly if realistic transit light curves are 
provided to support the claimed predictions. 

In the following, we assume that the early absorption fea- 
ture observed for HD 189733b could be caused by its magneto- 
sphere. To test this scenario, we need to estimate the stellar wind 
conditions at the planet's orbit at 0.031 AU as well as the spa- 
tial orientation of the nose of the magnetosphere relative to the 
planet- star line, taking into account the orbital velocity of the 
planet (IVidotto et al. ll2010h . HD 189733 is a main-sequence K 
star for which a coronal stellar wind is expected as revealed by 
the intense X-ray emission observed (ISanz-Forcada et aLll20 1 U 
iLecavelier et al. Il2012l) . To model the stellar wind formation, we 
used in a first step the Parker mode l which represents a good ap - 
proximation despite its limitations dParker II 19581 ICohen 11201 lb . 
X-ray observations show that the HD 189733 corona has a tem- 
perature in the range of a few 10 s - 10 6 K and a density up 
to ~ 10 10 crrT 3 (Sanz-Forcad a et al. 11201 lb . We accordingly as- 
sumed these properties as boundary conditions for the base of 
the stellar wind in the Parker model. Combining the isothermal 
Parker model with a hydrostatic distribution of the gas provides 
the stellar wind parameters (speed, density, and temperature) at 
the planet's orbit as summarized in Table 4. For the interplan- 
etary magnetic field (IMF) at the orbit of the exoplanet, we as- 
sumed a field strength in the range ~ 4 - 23 mG reported for 
the average field f rom polarimetry observations of HD 189733 
dFares et al. I l2010b . At the level of our exploratory work, we 
can free ourselves from accuracy concerns of the Parker model 
and accept any solution for the stellar wind properties within 
the range displayed in Table 4. Taking into account these wind 
parameters, the orbital motion of the planet, and the stellar coro- 
tation effect, we obtain a geometrical orientation of the magne- 
tosphere with a nose de flection angle in the range 10 - 30° from 
the exoplanet-star line ( Vidotto et al. 1 1201 ll) . 

To quantitatively estimate the opacity impact of the magne- 
tosphere on the stellar flux around transit, we first need to model 
the HD 189733b magnetosphere using the stellar parameters in 
the range described in Table 4. Here, we assume that the plane- 
tary magnetic field is a dipole, whose strength is a fre e parame- 
ter in the range < 14G dReiners & Christensen l uOlO). Because 
our goal is to derive the large-scale structure of the exoplanet's 
magnetosphere, we used a pa rticle-in-cell (PIC) 3D elect r omag- 
netic and rela tivistic code (Bara ka & Ben-Jaffel 1 120071 12011; 
IWodnickall2009h . which we describe in the following section. 



4.1 .1 . PIC simulation of magnetosphere and comparison 
with observations 

We employed a PIC cod e, originally developed and validated for 
Earth's magnetosphere (Buneman1 ll993t iBaraka & Ben-JaffeTI 



120071 120111) . which we extended to HD 189733b using the 
fact that the expec ted magnetospheric cavity is relatively small 
dWodnicka 1 120091) . Indeed, assuming a typical magnetic field 
strength of ~ 7G dReiners & Christensen I I20T0I) and a stellar 
wind of ~ 500km/s, the magnetopause position is estimated to 
be ~ 6 - 8 R p comparable to ~ 10 R p on Earth. This allowed us 
to use a relatively dense grid to describe the whole system with- 
out requiring very expensive computer resources. Electrons and 
ions are represented as macro-particles that contain a large num- 
ber of real particles. The PIC code solves Maxwell equations on 
the assumed grid 



dB 

— = - v xE, 
dt 



dE , 
e — = [i l V xB - J, 



(1) 
(2) 



where J — E (n,g,v; - n e q e v e ) is the current vector, and follows 
each macro-particle in the simulation box using the Newton- 
Lorentz motion equation, 



d{jmv) 
dt 



q(E + vxB), 



(3) 



accounts for the relativistic motion of 



where y = ^jj^Ty 
particles. The PIC code units are such that the dielectric 
constant e a = 1 and the magnetic permeability of vacuum 
Ho = 4. Because charge is conserved and magnetic fields are 
divergence-free, both Gauss laws appear as initial conditions. 
In other terms, since the Gauss laws are fulfilled at t=0, they 
will remain satisfied at any future time step o f the simulation 
dVillasenor & Buneman uneman 1 119931) . 

To avoid numerical i nstability, the Couran t condition cAt < 
should be fulfilled dCourant et al. I Il928l) : here the speed 

of light speed was taken to be c = 0.5, At is the step time, 
and Ar is the unform grid size. Plasma instabilities can be ef- 
ficiently reduced with t he strong condition on the plasma fre- 
quency Up. At < 0.25 dTskhakava et al. 1120071) . In addition, to 
avoid the classical problems of grid heating and instabilities, we 
ensured that the Debey length did not reach values below a crit- 
ical s ize defined by Ad > ^ dCai et al. Il200l iTskhakava et al. I 
2007). To enforce the shielding effect of charges on the Debye 
volume, we ass umed an initial particle de nsity of five pairs per 
simulation cell (Reimann & Faiansl 120021) . In the PIC code, the 
strength of the magnetic field of the planet is selected to ob- 
tain a magnetopause position at the desired value. An initial esti- 
mate was obtai ned from emp i rical models of the magnet opause 
position dVidotto et al. Il2010t iKhodachenko et al. II20T2I) . Non- 
reflecting boundary condi tions were applie d to the fields on all 
external facets of the box ( Lindm anll 1 9751) . Corotation with the 
planetary magnetic field was included by adding a corotational 
electric field as a boundary condition at the surface of the planet. 
Particles traveling down toward the surface are lost and thus re- 
moved from the simulation box. This explains the empty volume 
obtained around the planet's position (see Fig. 7). 

The simulation box has dimensions of 305 x 195 x 195 Ar 5 
along the OX-OY-OZ axes, each pixel having a uniform size in 
the range Ar = 0.2 - 0.5 R p depending on the size of the mag- 
netospheric cavity and limitations related to our computer re- 
sources. Simulations started with a total number of 5.5 x 10 7 
pairs of macro-particles in the box. The OX axis corresponds to 
the impinging stellar wind direction relative to the planet and 
makes the angle 8 with the planet-star line (e.g., Table 4). The 
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Table 4. HD 189733 stellar wind and magnetosphere parameters 

Density (cm -3 ) Temperature (K) Speed (km/s) Bimf (mG)'"' Magnetopause nose 8 angle' 
< 5 x 10 7 < 4 x 10 6 ~ 200 - 900 4-23 ~ (10 - 30)° 



Notes. All quantities are estimated at the exoplanet's orbital distance of 0.031 AU. The wind parameters are based on the Parker model and a 
hydrostatic distribution of isothermal gas density. The bottom conditions of the stellar c orona are defined from X-rays observations (see text), 
(a) The interplanetary magnetic field is based on polarimetry observation of HD 189733 dFares et aiT i 2010). (b) The angle is derived using an 
orbital velocity of ~ 152.5 km/s and a stellar corotation speed of ~ 28.5 km/s of the stellar plas ma and magnetic field at the distance of the planet 
dVidotto et al. ll201lh . A spin period of 1 1.9 days is assumed for the star dHenrv & Winn 1120081) . 




Fig. 7. 3D PIC code simulation of the HD 189733b magneto- 
sphere as seen from the bottom of the simulation box. The equa- 
torial plane is shown. We can clearly see the exoplanet's bow- 
shock and atmospheric nebula. Streamlines show the large-scale 
flow pattern, particulaly the corotation effect of the planet. This 
simulation was obtained for a stellar wind that has a sonic Mach 
number M s ~ 4.8, an Alfven Mach number Ma ~ 5.5, and a 
magnetosonic Mach number M m s ~ 3.6. The exoplanet is lo- 
cated at coordinates x = 120Ar, y = 97Ar, and z = 98Ar in the 
simulation box (reference frame defined by the star-planet line, 
the north-south magnetic poles, and the dawn-dusk directions). 
The intrinsic magnetic field of the planet is selected to obtain a 
stand-off distance of ~ 16.7^, which allows a good fit to the 
CII light curve (see Fig. 8 ). The grid has Ar units. The star is 
located toward the top of the page. 



OZ axis coincides with both the planet and star spin axes. The 
XY plane is taken as the equatorial plane of the planet that here 
concides with the orbital plane. Indeed, the HD 189733 spin 
axis is believed to be aligned with orbital plane of the planet 
(ICohen et al. 11201 ll) . The ion-to-electron mass ratio used in the 
code mi I me - 32 is high enough to obtain a good separation 
between opposite charges and helps complete the simulation i n 
a reasonable amount of steps (Gargat Tand Spitkovskvll2012l) . 
Finally, the code parameters selected thus far yield an ion skin 
depth A, ~ 12 Ar, a value that ensures that the mag netospheric 
cavity is properly sketched with the selected grid (Omid iet al. I 
20Q4: lMoritaka et al. 120121) . We stress here that our goal was not 
to resolve the detailed kinetic properties of any specific region of 



the magnetosphere but rather to model its large-scale geometry 
in a coherent way, which should help test its level of absorption 
during transit at specific ion resonance lines. 

We modeled the extended atmosphere as a cloud of pairs of 
macro-ion and macro-electron pairs that are randomly and con- 
tinuously injected to produce a uniform spherical flow around 
the planet that is confined in space and that has the same ki- 
netic temperature ~ 10 4 K, d ensity of few 10 7 cirT 3 , and radial 
speed ~ 5 km/s (lGuoll20TTl) . The spatial extent and the total 
content of the resulting nebula were free parameters that were 
to be fixed by the transit absorption depth. Ideally, one should 
include an altitude-dependent ionospheric distri bution, but that 
is bey ond the scope of this preliminary study dJohansson et al. I 
2009). Initial conditions require Maxwell distributions for the 
velocity field of both macro-ions and macro-electrons at their 
corresponding temperatures. Stellar wind was injected versus 
time from the YZ plane located on the top edge of the box in the 
OX axis (see Figure 7). The planet's dipole magnetic field was 
switched on smoothly until it reached its maximum value. The 
system was then left evolving to reach a steady state in which the 
observable large-scale structures were not changing anymore. 
Post-processing of simulation results allowed us to derive all 
moments of the macro-particles velocity distributions anywhere 
in the box. Plasma number density, temperature, and bulk veloc- 
ity distributions were thus calculated, allowing us estimate the 
gas opacity of the magnetosphere at the desired spectral line and 
its obscuration of the stellar disk flux as the planet moves along 
its orbit from ingress to egress. In Figure 7, a typical magneto- 
spheric configuration is shown in the equatorial plane. We can 
easily distinguish the planetary atmospheric nebula, the plasma 
equatorial sheet, and the upstream bow shock. Streamlines of the 
bulk speed in that plane show the corotation effect on the plasma 
located near the planet. Here, it is important to stress that the 
equatorial plasma and current sheet (also called magnetodisk) is 
naturally formed with the PIC code, a consistency that cannot be 
met in MHD simulations. 

In a first step of the sim ulations, assuming a surface mag- 
netic field strength of ~ 7 G dReiners & Christensen 1 120101) for 
the planet and using the fast stellar wind solution as derived from 
X-ray observations (~ 500 km/s ), a magnetospheric configura- 
tion with a nose angle of 10 - 30° from the star-planet line pro- 
duces only a poor fit to the C II light curve (assuming C/H solar 
abundance; see Fig. 8). It is important to notice that we do obtain 
an asymmetric transit light curve and early-ingress absorption, 
but later than observed. The model misfit tends to show that the 
magnetopause (MP) position obtained by the model (~ 7 R p ) is 
too close to the planet to enable a fit to the observed C II early- 
ingress absorption seen in HST orbit 1 (see Fig. 8). 

To improve our diagnostic, we therefore focused on the mag- 
netopause position that allowed to fit the observed early-ingress 
absorption. In practice, we varied the ram pressure of the stel- 
lar wind until a good fit was obtained for the C II light curve, 
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HD 189733b Transit (C II 133.5 nm) 
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Fig. 8. 3D PIC simulation Fit to HD 189733b C II 133.5 nm light 
curve. The simulations are obtained for the stellar plasma pa- 
ramete rs indicated in Table 4 . A solar C/H= 2.45 x 10 is as- 
sumed ( Asplu nd et al. 1 12005). (Dashed) Model assuming a mag- 
netic field strength of 7 G for the exoplanet and a fast stellar 
wind (speed ~ 500km/s, density ~ 6 x 10 6 ctrT 3 , and temper- 
ature ~ 4 x 10 5 K). (Solid) Using a magnetic field strength of 
~ 5.3 G for HD 189733b, and stellar wind speed of ~ 250km/s, 
temperature ~ 1 .2 x 10 5 K, and density ~ 6.3 x 10 6 cm -3 lead to a 
better fit corresponding to a MP position of ~ 16.7 upstream. 
This solution is not unique (see text). 



keeping the magnetic field strength of the exoplanet assumed in 
the PIC code unchanged. As shown in Figure 8, a good fit is ob- 
tained with an MP upstream stand-off distance of ~ 16.7 R ; ,. It is 
interesting to notice that in addition to the transit signature due 
to the planetary atmosphere (at ±3000 s from mid-transit), early 
absorption by stellar plasma in the bow-shock region clearly ap- 
pears on ingress (around 10,000 s before mid-transit), while a 
weak absorption also appears on late egress (~4000- 10,000 s af- 
ter mid-transit), probably from plasma in the magnetotail (see 
Fig. 8). We emphasize that because the magnetopause is a stag- 
nation region that results from total pressure balance between 
the star and the planet, other solutions for the planetary magnetic 
field and stellar wind parameters may lead to the same MP po- 
sition. Furthermore, because the stellar parameters are not accu- 
rately known, the orientation of the magnetosphere is not unique, 
which leads to a degeneracy in the possible solutions of the pro- 
jected size of the magnetospheric cavity. 

For example, for a stellar wind speed of ~ 250km/s, a tem- 
perature of ~ 1 .2 x 10 5 K, and a density ~ 6.3 x 10 6 cirT 3 at the 
planet orbit, the planet's magnetic field should be ~ 5.3 G to fit 
the MP position at ~ 16.7 R p . Assuming the same stellar wind 
parameters, but with a speed of ~ 450km/s, our model simula- 
tion requires an exoplanet magnetic field of ~ 7.2 G to obtain a 
good fit to the C II transit light curve. As a final example of a 
good fit, stellar wind parameters with a speed of ~ 250 km/s, a 
temperature of ~ 5.0 x 10 5 K, and a density ~ 1.3 x 10 7 crrT 3 
at the planet orbit along with an exoplanet magnetic field of 
~ 12.5 G also give a satisfactory fit to the C II transit light curve. 



All these results can be easily understood from the total pres- 
sure balance that occurs at the MP. Indeed, if the stellar wind 
total pressure is changed for the same MP position, the total 
planetary pressure should be modified on the opposite side to 
maintain the pressure balance. In that way, an excessively strong 
planetary magnetic field could be required to fit the MP posi- 
tion at ~ 16.7 R p for some stellar wind parameters. For ex- 
ample, for a stellar wind temperature of ~ 10 6 K, a density of 
~ 1.0 x 10 7 crrT 3 , and a speed of ~ 250 km/s, a good fit could be 
obtained but with a planet magnetic field of ~ 16.2 G. Finally, 
we emphasize that we could find models that produce a magne- 
tosphere with the required MP position, but that do not fit the 
C II transit light curve. A full sensitivity study versus the stellar 
wind and planetary atmosphere parameters is needed, but this is 
beyond the scope of the present study. 

In practice, to improve the present diagnostic of a magneto- 
spheric signature during transit, we need to obtain a better cover- 
age of the transit light curve of strong resonance lines of several 
ionized species, particularly on ingress and egress (see Fig. 8). 
Additionally, to account for potential orbital phase variations in 
the stellar wind, we may need models that include the full stellar 
corona-exoplanet system in the simulation b ox to better uncove r 
any star-planet electromagnetic coupling (ICohen et al. 1 1201 Tb . 
For all these reasons, only a sensitivity study may help separate 
the different effects and properly constrain the exoplanet mag- 
netic field if a full transit light curve is confirmed with the ex- 
pected strong early-ingress or weak late-egress absorptions. To 
that end, recent studies on the solar wind interaction with the 
magnetized local interstellar medium can provide a good refer- 
ence for a future approach to uncove r the plasma environment 
and m agnetic field of HD 189733b dBen-Jaffel & Ratkiewiczl 
[2012h . 



5. Summary and conclusions 

Using archival HST/COS G130M observations, we have posi- 
tively detected neutral oxygen ~ 7% absorption during the tran- 
sit of exoplanet HD 189733b. We also reported a tentative de- 
tection of singly ionized carbon. If sketched as a compact disk, 
this detection of OI reveals a distribution of atoms extending 
to about ~ 1.1 R p inside a Roche lobe radius of ~ 2.84 R,,. 
Assuming a mean temperature of ~ (8 - 1 2) x 1 3 K for the upper 
atmosph ere of the exoplanet, we used the hydrodynamic model 
of lGuol d201 lb to derive that for solar abundances with an OI 
density column of ~ 8 x 10 15 cm -2 , an attenuation of ~ 3.5% is 
obtained for the OI 130.6 nm line during transit, at least ~ 2.5<x 
below the observed absorption. Invoking up to 40 times solar 
abundances with an OI density column of ~ 3.3 x 10 17 cm" 2 , 
produces an attenuation of ~ 5.3% during transit that is closer to 
observations. These very high density columns of heavy atoms 
are difficult to justify because they would require substantial ad- 
ditional energy in the system to compensate for the reduced scale 
height. In contrast, we show that including a super-thermal O I 
layer embedded on top of the atmosphere reduces the need for 
highly super-solar abundances. For example, with twice solar 
abundaces and a hot O I layer confined above ~ 1.5 R p with 
an effective temperature of Toi ~ 8.4 x 10 4 K, a transit ab- 
sorption of ~ 4.3% is produced that is only about ~ Icr be- 
low the observed value. Our preliminary results tend to con- 
firm the general conclusion that similar to HD 209458b, super- 
solar abundances with thermal broadening as well as super- 
thermal broadening probably affect th e FUV transit of hot ex - 
oplanets dBen-Jaffel & Hosseini Il2010t iKoskinen et al. Il2012bh . 
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High-resolution line profiles observations of key FUV lines must 
be obtained with high S/N to separate the different effects. The 
O I detection reported here opens up new perspectives to use the 
high sensitivity of the HST/COS instrument to detect oxygen 
in extrasolar systems in which the stellar O I lines are thin and 
bright enough compared to the broad geocoronal line that results 
from the extended sky background filling the 2.5" aperture. 

Furthermore, we found a peculiar signature in the CII 
133.5 nm transit absorption by HD 189733b that shows differ- 
ent properties and a dependence on orbital phase, which appears 
to be as an early-ingress absorption. Combined multi-spectral 
analysis and a comparison of line profiles support the reality of 
the feature. However, because the star is relatively active (e.g., 
Fig. 6), both the transit and the early-ingress detections require 
confirmation in the future to determine any stellar variability that 
cannot be controlled from the present analysis of only one single 
transit observation. 

Assuming the extra absorption is real, we used the Parker 
model for the stellar wind and a PIC code simulation of the 
interaction of the magnetized and extended thermosphere of 
HD 189733b with the impinging plasma wind. First, our analysis 
showed that a magnetosphere forms with a bow shock and a nose 
that is oriented almost ~ 10 - 30° from the star-planet line for 
most cases of stellar wind parameters. Second, the simulations 
revealed that a magnetopause stand-off distance of ~ 16.7 R p 
is required to obtain a satisfactory fit to the transit light curve 
observed for C II (e.g. Fig. 8). However, our preliminary assess- 
ment of the PIC simulation results showed that many solutions 
exist for the stellar wind parameters and the magnetic field of 
the exoplanet that fit the C II transit curve. Several reasons could 
explain the origin of this degeneracy, namely the limited cover- 
age of the transit event by the HST/COS observations, the sig- 
nal variability (stellar or instrumental), and finally the nature of 
magnetospheric cavity. Indeed, because the magnetopause is a 
stagnation region where a total pressure balance occurs between 
stellar and planetary plasmas and fields, no unique solution ex- 
ists for either the stellar wind parameters or the magnetic field 
strength of the exoplanet that provide a prescribed magnetopause 
stand-off distance. For example, we are able to find a good fit for 
the CII transit light curve with a speed of 250km/s, a tempera- 
ture of ~ 1 .2 x 10 5 K, and a density ~ 6.3 x 10 6 cuT 3 at the planet 
orbit for the stellar wind, and a magnetic field strength of ~ 5.3 G 
for the exoplanet magnetic field. However, other solutions exist 
(see previous section). 

In all cases, whether the C II peculiar feature is real or not, 
our PIC simulations showed that the C II spatial distribution that 
contributes to the observed transit absorption is not solely re- 
lated to the planet's extended atmosphere but rather is dependent 
on the plasma configuration that results from the complex inter- 
action between the stellar wind, the exoplanet atmosphere, and 
the local planetary and stellar magnetic fields. New HST FUV 
observations are urgently needed to confirm our results and ex- 
tend the coverage of the early-ingress and late-egress behavior in 
the C II 133.5 nm line absorption. Our quick comparison of PIC 
simulations of HD 189733b's magnetosphere with transit obser- 
vations showed that fundamental properties of the plasma con- 
figuration around the exoplanet may be accurately constrained 
through the specific shape of the light curve versus wavelength 
bands either from a single line or using multiple lines. 
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